syms x y q12_re  q12_im

lambda = 1+2i;

w = i;
w = 1;

r = solve(real(lambda^2) + sqrt(q12_re^2+q12_im^2) + real(w*(q12_re+i*q12_im)) -x^2+y^2,...
    abs(lambda)^4 +2*abs(lambda)^2*sqrt(q12_re^2+q12_im^2) - (x^2+y^2)^2 + ...
    real(lambda)^2*real(w)*q12_re - real(lambda)^2*imag(w)*q12_im - 2*real(lambda)*imag(lambda)*real(w)*q12_im - 2*real(lambda)*imag(lambda)*imag(w)*q12_re - imag(lambda)^2*real(w)*q12_re+imag(lambda)^2*imag(w)*q12_im     ...
 ,...
    q12_re, q12_im);
